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Abstract 

Short-time dynamics technique is used to study the relaxation process for the (2+1)- 
dimensional critical gluodynamics of the SU(2) lattice gauge theory. A generalized 
class of heat-bath-inspired updating algorithms was employed during the short-time 
regime of the dynamic evolution for performance comparison. The static and dy- 
namic critical exponents of the theory were measured, serving as a dynamic bench- 
mark for algorithmic efficiency Our results are in agreement with predictions from 
universality hypothesis and suggest that there is an underlying universal dynamics 
shared by the analyzed algorithms. 
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1 Introduction 

The pioneering study on equilibrium critical phenomena addressed by Fisher 
et. al. through finite-size-scaling relations [1] was soon extended to include a 
description of dynamic relaxation processes 02] ■ Thus, it was demonstrated 
that in many systems an universal behavior holds even far from equilibrium 
|4T5ll6] . In lattice-gauge theories, Monte Carlo relaxation dynamics was investi- 
gated by Okano et. al. in seminal works [7f8|9yi0] where critical-exponents were 
computed. These studies reinforced the hypothesis that universality among 
spin-systems and gauge-theories jTT] is dynamically realized. In addition, pre- 
dictions from universality were checked in lattice-gauge theories beyond the 
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determination of critical-exponents. For instance, it was confirmed that, near 
the deconfinement phase-transition, the dynamically generated spectrum of 
gluonic-screening masses obeys universal ratios [T2lll3|ll4] . 

Equilibrium Monte Carlo simulations allow for non-perturbative numerical 
calculations, so they are essential for the study of some fundamental theoret- 
ical issues on phase-transitions of lattice-gauge theories |15j . However, when 
approaching the thermodynamical limit, those simulations have their thermal- 
ization efficiency increasingly afflicted by the so called critical-slowing-down 
effect |16j . which makes it costly to compute independent gauge configura- 
tions. A possible way to circunvent that effect comes from Short-time dy- 
namic techniques |8fTU] . which enables critical gluodynamics to be efficiently 
exploited by local [T7|18lll9j or global [20j updating algorithms. Therefore, to 
devise improved algorithms [21.22.23j by better understanding their equilibra- 
tion features — e.g. (over)relaxation dynamics — may be worthy for many 
applications. 

In this context, we extend our previous comparative study over a class of 
improved heat-bath-inspired thermalization algorithms [221123] to the non- 
equilibrium regime, by using a (2+l)-dimensional SU(2) lattice-gauge theory 
at critical temperature. It is known that dynamic-relaxation exponents 9 and 
z — namely in contradistinction to static ones: a, f3, 7, 5, v — are generally 
dependent on the very dynamics of thermalization; so, they would serve as 
discriminants for algorithmic classes. The article is organized as follows: in 
Section 2, general properties of short-time critical dynamics are reviewed; an 
overview of usual heat-bath algorithm and its generalizations is the theme of 
Section 3; the setup of our simulations and the description of a new procedure 
to implement sharp initial states on lattice-gauge theories is given in Section 
4; general data analysis and simulation results are provided on Section 5; we 
summarize our findings in Section 6. 



2 Short-time critical dynamics 

Renormalization group techniques allowed Janssen et al |3] to show that un- 
der suitable initial conditions some magnetic systems, after being suddenly 
quenched to critical-temperature T c , present universal-scaling behavior even 
at early evolution. This phenomenon is observed just after a microscopic tran- 
sient time-scale has elapsed and lasts for a macroscopic relaxation-period 
tmacro before thermalization is established. 

During the short-time process of relaxation, when the system is driven to 
equilibrium by a Monte Carlo stochastic-dynamics, just mild finite-size-effects 
[7] are noticeable. Then, reliable information can be extracted from the re- 
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suiting time-series by considering some generalized-scaling relations for usual 
observables. For instance [7f5] , the kth moment of magnetization obeys 



M w (t, r, L, m ) = b~ k ^ u M^ (jr g t, b 1/u r, b^L, b x °m ) , (1) 

where itlq is the initial magnetization, b is an arbitrary spatial scaling factor, 
L is the system size, t is the time-evolution parameter, and the reduced- 
temperature is r = (T — T c ) /T c . While /3 and v are the well-known critical- 
exponents of equilibrium (i.e. static), z and xq are relaxation-dependent (i.e 
dynamical) exponents, and both groups label universality classes. 

Considering a hot initial-state with sharply-defined (and small) magnetization 
mo, it was shown |4"|5] that for a quenched system, some simpler power- law 
relations hold for the magnetization 

M(t,m ) «m t fl , (2) 
and its higher-moment 

M (2) (t) _ t (dr-2P/u)/Z } (3) 



where 9 = (xq — /3/u) / z and d denotes the spatial dimensionality. Analo- 
gously, the temporal-autocorrelation A (t) for the magnetization evolves as 

A(t)~t e - d/z . (4) 



Consequently, the critical phenomenology of a system can be studied by ob- 
serving the temporal-scaling behavior — up to the time-scale t ~ mQ 2 ' 10 
[I] - - of some of its observables, during the early relaxation (i.e. short-time 
regime) of the dynamic evolution. This technique allows for the determi- 
nation of equilibrium and non-equilibrium critical-exponents [3j without se- 
vere critical-slowing-down restrictions [7]. Hence, averages in the dynamic 
approach are taken over independent samples, without relying on ergodical 
time- averaged measurements [8]. 



3 Heat-Bath inspired algorithms 



Thermalization is the process responsible for producing independent config- 
urations during Monte Carlo simulations. When the quenched approxima- 
tion is used in gauge theories, by setting the fermion determinant to unit, it 
is possible to apply local algorithms such as Metropolis or heat-bath (HB) 
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[TTfT8] as efficient first-choices for thermalization. However, when a critical- 
point is approached, as in the continuum limit of the theory, simulations suffer 
from critical-slowing-down phenomenon [24]. This drastically increases corre- 
lations among successive field-configurations, which produces integrated auto- 
correlation times Ti n t raising as a power of the lattice side L. Thereby, it 
induces Monte Carlo statistical errors to diverge as O (v / 2t^^ |16| . 

In order to circumvent the critical-slowing-down, the standard HB-algorithm 
and micro-canonical updates are combined, which allows for improved gen- 
eration of independent samples |24j . In particular, in the hybrid version of 
overrelaxed algorithms, m energy-conserving microcanonical update sweeps 
are done after a standard local ergodic update. 

For a SU(2) lattice-gauge theory, the action can be factorized as a sum of 
single-link actions Si-u n k, which may be written down as 

Si-imk = —^Tr[Up(x)Hp(x)], (5) 

where U ^ (x) G SU (2) is a gauge-link, and H n (x) — the effective-magnetic 
field, written as a sum over staples — is proportional to an SU(2) matrix. As 
a useful notation consider H n (x) = (x) H a (x) with (x) G SU (2) and 

(x) = ^JdetH u (x). 

Then, by using Eq.(|5]) and the invariance of the group measure under group 
multiplication, the HB update is obtained 

where V = VqI + i-v ■ a G SU (2) is randomly generated by choosing Vq 
according to the distribution y/l — f 2 exp ((3Nvq) dvo and v pointing along a 
uniformly chosen random direction in IR 3 , under the constraint | V" | 2 = v g + 
v-v = 1. 

We have proposed in [22.23J a modified HB algorithm (MHB) in which the 
generation of the updating matrix V is carried out as usual, except for the 
additional step 

• Transform the new vector- components of V as v — >■ —sgn(v ■ w)v, 

where W = WqI + i ■ w ■ a = U° ld (x) H u (x), and sgn is the sign function. 

This may be also thought as an ergodical modification of the overheat-bath 
(OH), an algorithm devised some years ago [21] • The underlying idea in both 
cases is to incorporate a micro-canonical move [23] into the heat-bath step. 
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The difference is that while in MHB the vector v is randomly set — except 
for its sign, which is determined according to the aforementioned rule — in 
the OH algorithm, one deterministically sets v = —w — without obeying a 



uniform distribution for v — and renormalizes it as \v\ — y 1 — Vq. Therefore, 
OH incorporates a micro-canonical move in an exact algorithm, but it would 
not be ergodic 



Despite being slightly more computer-consuming than OH, our modification 
implements a micro-canonical move that explicitly preserves ergodicity and 
allows for a reduction of about 20% in statistical errors, at the same compu- 
tational cost, when compared with the usual HB algorithm |22|23] . 



4 Simulation setup 



Lattice gauge theories in d-dimensions can be easily turned to finite-temperature 
formalism [8|12lll3f25| by constraining the euclidean spacetime volume to V = 
L d ~ 1 L t — under the assumption L t <^ L s — where L s and L t are the spatial 
and temporal lattice sides. The equilibrium temperature is given by T _1 = 
a ■ L t , where a is the physical lattice spacing. 

The Polyakov loop P$ on site x is a useful quantity whose spatially averaged 
trace J2gTr(P%) is the order parameter of the deconfinement phase-transition. 
It is defined as an ordered product of gauge-links Uq (xq,x) in the temporal 
direction. Here we use, for the particular SU (2) case, a non-equilibrium time- 
dependent definition 

L 3 [t] = Tr {P s } [t] = Tr { \ f[ U {x , x) ) [t] , (7) 

1 1 zo=o J 



where [t] denotes a given instant in the Monte Carlo time. 

A strict analogy with the magnetization in Eq.(|2]) is addressed to gauge the- 
ories |8j by defining a time-dependent effective magnetization M as 

^(*) = (^EM«l) - (8) 

' s x ' sample 



where (• • •) sample stands for an average taken over simultaneous samples, at an 
instant t. Similarly, Eq.([3]) is associated to the second-moment of M by the 
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relation 



M<® (t) = 



sample 



(9) 



and Eq.Q is related to 



A(t) = 




sample 



(10) 



For the purpose of performing short-time dynamic simulations, it is necessary 
to prepare hot-initial states (T = oo) with a sharply-defined (small) effective 
magnetization (M m tuq). This may be achieved by the procedure 



Setting a sharp initial-state 

(1) Restart the random-number generator using a new seed. 

(2) Set all components of each SU (2) gauge-link randomly. 

(3) For a maximum number N steps of loops do 

(4) Evaluate M (0) as in Eq.([8]). If \\M (0) — rrio\\ < b~m to i erance start step-8, 
else do 

(5) Choose up to N S i tes lattice-sites at random positions (xq,x). Then, update 
their temporal gauge-links by Uq (xo, x) — > —sgn [M (0) — mo] PgUo (xq, x). 
Where P% is as in Eq.(j7]). 

(6) Return to step-4. 

(7) Return to step-3. 

(8) Perform thermalization. 

In our simulation-setup, parameters were tuned to N steps ~ O (10 4 ) and 
N S i tes ~ 0.2 x Lg' 1 , which allowed us to shortly obtain <5m to / erance /mo ~ 
O (10~ 3 ) , without introducing any spurious spatiotemporal correlations among 
samples. 



5 Numerical results 

We have used lattice sizes up to 128 2 x 2 in our simulations. Initial lattice 
configurations were set using the previously described algorithm. Despite the 
dynamic critical exponents in Eq.(|2]), Eq.Q, and Eq.Q being rigorously de- 
fined just for m = 0, the minimal initial magnetization we considered was set 
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to be mo = [400(1)] • 10 5 — i.e., no extrapolation to m — > was attempted 
— which allowed us to keep better signal-to-noise ratios. 



To thermalize the gauge-fields, we employed the aforementioned heat-bath-like 
algorithms and the standard Wilson action at critical-coupling = 3.4505 
|10j . We have prepared 50000 independent initial samples for simulations, 
whose measured observables [Eqs. ^ — ( 
estimate errors — for each lattice sweep. The Monte Carlo temporal evolution 
was followed up to 500 steps, as seen in Figure(l): Left panels. 



10)] were grouped in 10 blocks — to 



The best-fit range for each observable was determinated by searching for the 
largest plateaus in the time-interval (^,500], within the minimal x 2 /dof, 
which were produced as outcomes for previous power-law fits, Figure(l): Right 
panels. Within the established stable-regions we accurately fit the correspond- 
ing critical exponents. This methodology also allowed us to unveil the ther- 
malization efficiency of the employed algorithms by their direct comparison 
during relaxation. 



It is interesting to note that average numerical-values from data of time-series 
produced by MHB-evolution are nicely located between the ones obtained by 
HB and OH algorithms, Figure(l): Left panels. Also, despite their different 
inner-dynamics, the algorithms we implemented reproduced self-consistently 
- up to good numerical approximation — the same set of static and dynamic 
critical-exponents, as seen in Table(l). 



The numerical results for critical exponents agree with previous studies |8f9fT0f2 
that investigated the universality-hypothesis |TT] between SU(2) lattice gauge 
theory and the 2d-Ising model. Also, upon closer examination for similarities 
among the time-series produced by different algorithms — and their dynamic 
9 and z exponents — there is a strong indication that all heat-bath-inspired 
algorithms analyzed here share the same universal non-equilibrium-relaxation 
dynamics/ 1 "] 



This may be a hint that effects of some hypothetical ergodicity violations, pre- 
viously conjectured in |23j . would be less harmful than what is expected for some 
simulations [12] using the OH algorithm. 
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Figure 1. Monte Carlo time-evolution of observables in short-time relaxation dynam- 
ics for different algorithms. (Panel A) and (Panel B) are plots for the first and 
second moments of the magnetization, (Panel C) is the temporal auto-correlation 
of the magnetization. Error bars are simple Monte Carlo standard deviations, they 
have been omitted (on Panel A) for better visualization. 

Detection of the fitting-range [t^,t] by analysis of the x 2 /dof vs. t produced as an 
outcome from power-law fit anzates for different observables. (Panel D) and (Panel 
E) are plots for the fit-quality for the first and second moments of the magnetization, 
(Panel F) is the fit-quality for the temporal auto-correlation of the magnetization. 
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Algorithm 


Q 
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p/u 


UD 

o 1 3 
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U. llU^l j 


MHB 


0.1843(7) 


2.141(4) 


0.114(5) 


OH 


0.1955(1) 


2.094(7) 


0.132(7) 


Ising-2d (HB) 


0.191(1) 


2.155(3) 


1/8 



Table 1 

Static and dynamic critical-exponents obtained with each algorithm employed for 
time-evolution in our Monte Carlo simulations. The last line presents the same 
quantities obtained for the two-dimensional Ising model [10J. 



6 Summary 



Our results extend previous studies on short-time dynamics of lattice gauge 
theories [8, 9,10,20] to a whole class of heat-bath-inspired thermalization algo- 
rithms. This would be seem as a self-consistent cross-check of the universality 
hypothesis |TT] holding far-from-equilibrium. Here, this is verified for a set of 
thermalization algorithms that embraces different relaxation dynamics. 

Notwithstanding the particularities of each overrelaxation-dynamics embed- 
ded in the algorithms we analyzed, they seem to display an underlying uni- 
versal dynamics, a fact that is also corroborated by the numerical values ob- 
tained for the dynamical critical-exponents. Thus, in this comparative study 
for benchmarking thermalization algorithms, when focusing on finding a bet- 
ter balance among most desirable features for updating algorithms — e.g., 
faster decorrelation of samples and higher signal-to-noise levels at same com- 
puter cost — it is observed that MHB presents the most desirable algorithmic 
realization. 

The methods employed in this article can also be easily adapted to dynamic 
analysis of gauge theories for larger unitary, symplectic, or exceptional gauge 
groups [12] or even to O (N) continuous-spin systems [22] ■ That may bring 
some new perspectives for non-equilibrium comparative studies on universal- 
ity in full-QCD [25]. Also, some general constraints for the infrared gluon- 
propagator — written in close similarity with magnetization-moments [26] - 
would be investigated by short-time dynamic simulations, though with ame- 
liorated finite-size and critical-slowing-down effects. 
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